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Abstract 

The improved two-scale model with symmetric Lund distribution function is 
^ ■ used to perform the fit to the semi-inclusive deep-inelastic scattering (SIDIS) data 

on nuclear targets of HERMES experiment at DESY. The ratio of hadron multi- 
plicity on nuclear target to the deuterium one is chosen as observable, as usually. 
^ Some peculiarities of the fit are discussed. The two-parameter's fit gives satisfactory 

agreement with the data. The results of the present and previous fits are compared. 
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1 Introduction 

Hadronic reactions in a nuclear medium 
can shed additional light on the hadroni- 
zation process. In comparison with other 
reactions electroproduction has the virtue 
that energy and momentum of the struck 
parton are well determined, as they are tag- 
ged by the scattered lepton. Study of had- 
ron production in SIDIS on nuclear targets 
offers an opportunity to investigate the qu- 
ark (string, color dipole) propagation in nu- 
clear matter and the space-time evolution 
of the hadronization process. For this pur- 
pose we investigate the nuclear attenuation 
(NA), which is a ratio of differential hadron 
multiplicity on a nucleus to that on deu- 
terium: 



possible to sum analytically over all ranks 
of hadrons carrying momentum z. Unfor- 
tunately, more realistic symmetric Lund dis- 
tribution function it is not possible to sum 
over all ranks analytically. For this reason 
it was not included in the fit. 

The aim of this work is a fitting similar 
to that which was done in Ref. [3] with us- 
ing the symmetric Lund distribution func- 
tion and a comparison of the results of fit 
with the results of the previous one. 

The paper is organized as follows. In 
section 2 some details of the fitting pro- 
cedure with symmetric Lund distribution 
function are discussed. In section 3 we com- 
pare results of the present fit with the pre- 
vious one and discuss them. Conclusions 
are given in section 4. 
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where z = E^/v, Eh and v are energies of 
the final hadron and virtual photon, N h (v, z) 
is the number of semi-inclusive hadrons at 
given v and z and N e (v) is the number 
of inclusive DIS leptons at given v. Sub- 
scripts A (D) denote that reaction takes 
place on nucleus (deuterium). In the above 
formula more variables like -Q 2 (the pho- 
ton virtuality) and p t (hadron transverse 
momentum) over which the NA is aver- 
aged, are not written. The simple version 
of the string model, so called two-scale mo- 
del (TSM), was introduced in Ref. pQ. In 
Ref. [2] improved version of TSM (ITSM) 
was proposed. In our previous work [3] 
ITSM was used for a fitting of the SIDIS 
data of HERMES experiment on nuclear 
targets [!4J. For the fit we used those ver- 
sions of the distribution function, which 
possible to present in the finite form. In re- 
sult, two versions were used: (i) the leading 
hadron distribution function, which corre- 
sponds to the term with rank equal one 
of the distribution function; and (ii) the 
standard Lund distribution function, which 



2 Some features of fit 

In the string models, for the construc- 
tion of fragmentation functions, the scaling 
functions f(z) are introduced. For calcula- 
tions we use [5J E] : 

(i) standard Lund scaling function 

f(z) = (l + C)(l-zf, 

where C is the parameter which controls 
the steepness of the standard Lund frag- 
mentation function (C = 0.3); 

(ii) symmetric Lund scaling function 

f(z) = Nz-\l - z) a exp(-bm 2 ± /z) } 

where a and b are parameters of the model 
(a = 0.3, b = 0.58GeV~ 2 ), m± = yjm? h + p\ 
is the transverse mass of final hadron, N 
is normalization factor. We begin by con- 
sidering the distribution of the constituent 
formation lengths I of hadrons carrying frac- 
tional energy z. The simplified version with- 
out identifying flavor of partons has the 
form: 

D h c (L,z,l) = (f(z)5(l-L + zL) + 

n 

Y,D>\(L,z,l))6{l)e(L-zL-l) , 

i=2 



where L = v/k is the full hadronization 
length, k is the string tension (re = 1 GeV/fm). 
The functions D^(L,z,l) are distributions 
of the constituent formation length / of the 
rank i hadrons carrying fractional energy z. 
For calculations we use the average value of 
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Figure 1: x 2 as a function of iV^. Left 
part presents results of fit with standard 
Lund distribution function and right part 
results of fit with symmetric Lund distri- 
bution function. Panels (a), (b); (c), (d); 
(e), (f) correspond NDF 1; NDF 2; NDF 
3, respectively. Filled circles, open circles, 
filled triangles, open triangles are results of 
fit with a str corresponging eqs.(7), (8), (9), 
(10) of Ref. [3j. 

the formation length t c =< I >: 

roo /*oo 

t c = ldlD h c {L,z,l)/ dlD*(L,z,l) . 
Jo Jo 

The same versions of string-nucleon cross 
sections, a str , and the same set of the nu- 
clear density functions (NDF) as in Ref. [3] 
are used for fit. The semi- inclusive data jl] 
of HERMES experiment on four nuclear 
targets (helium, neon, krypton, xenon) and 
deuterium are used. As in previous fit we 



use one and two dimensional data for pi- 
ons, all together 332 experimental points. 
The fit was performed to tune two para- 



a str \l) 


NDF 


a q mb 


c 


x 2 


1 


0.00 ±0.01 


0.365 ±0.013 


1.48 


2 


0.00 ±0.01 


0.336 ±0.011 


1.32 


3 


0.00 ±0.01 


0.307 ±0.012 


1.30 


a str (8) 


NDF 


a q mb 


c 


x 2 


1 


1.74 ±0.16 


0.172 ±0.014 


0.81 


2 


1.88 ±0.15 


0.163 ±0.011 


0.71 


3 


1.84 ±0.15 


0.140 ±0.009 


0.74 


a str (9) 


NDF 


o q mb 


c 


x 2 


1 


0.00 ±0.01 


0.112 ±0.009 


0.75 


2 


0.00 ±0.02 


0.088 ±0.009 


0.66 


3 


0.00 ±0.02 


0.059 ±0.009 


0.66 


a str (10) 


NDF 


a q mb 


c 


x 2 


1 


2.17 ±0.13 


0.098 ±0.010 


0.61 


2 


2.34 ±0.14 


0.088 ±0.010 


0.54 


3 


2.43 ±0.15 


0.069 ±0.010 


0.59 



Table 1: Values of fitting parameters and 
X 2 in case of Symmetric Lund distribu- 
tion function and total errors. Numbers 
in parentheses indicate the corresponding 
equations from the Ref. [3]. 



meters: the initial value of string-nucleon 
cross section a q and coefficient c. 

The quantitative criterium x 2 was used: 
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where n exp and n par are numbers of exper- 
imental points and parameters; R^(theor) 
is the theoretical value for ratio at given 



2 The one (two) dimensional data means that 
data are presented as functions of one (two) vari- 
ables. 



point; R\ 4 {exp) and AR^iexp) are exper- 
imental value of R\j and its uncertainty at 
given point. 

Let's consider the features of fitting with 
symmetric Lund scaling function. The cor- 
responding distribution function it is im- 
possible to analytically sum over all ranks 
of hadrons as in case of standard Lund scal- 
ing function, therefore we must restrict our- 
self by finite sum over ranks of hadrons. 
The recursion equation from Ref. [7] is used 
for calculation of the distribution functions. 
We performed fits with the sums of ranks 
until n = Nh a d and convinced, that begin- 
ning from Nh a( i = 5 the values of parame- 
ters and x 2 do not change essentially. The 
close result is obtained for the case of stan- 
dard Lund scaling function. In this case the 
results for finite sums for N^ad > 5 prac- 
tically coincide with the result for infinite 
sum. Results of fits are presented in Fig. 1. 
In the left part of Fig. 1 (panels (a), (c), 
(e)) the results of fit with standard Lund 
distribution function and in the right part 
of Fig. 1 (panels (b), (d), (f)) the results of 
fit with symmetric Lund distribution func- 
tion are presented. Panels (a), (b); (c), (d); 
(e), (f) correspond NDF 1; NDF 2; NDF 
3, respectively. We would like to remind 
that NDF 1; NDF 2; NDF 3, correspond 
eqs. (14); (15); (16) of Ref. 0. Filled 
circles, open circles, filled triangles, open 
triangles are results of fit with a str corre- 
sponging eqs. (7), (8), (9), (10) of Ref. [3]. 

3 Comparison with data 
and discussion 

The results of fit with symmetric Lund 
distribution function and iV/^ = 10 are 
presented in the Table [TJ Easily to see 
that they describe data on quantitative level 
Minimum \ 2 = 0.54 is obtained for <j str (10) 
and NDF 2 from Ref. [3] at values of pa- 
rameters equal o q = 2.34 ± 0.14 mb and 
c = 0.088 ±0.01. Comparison with the Ta- 
bles 1 and 2 from Ref. [3] shows that values 
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Figure 2: The two dimensional data. The 
ratio R 1 }^ for charged pions on A He (panels 
a, b) and 20 Ne (c, d) nuclei function 
of v (left panels) and z (right panels). Ex- 
perimental points from Refs. jl], [8]. 
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Figure 3: The same as described in the cap- 
tion of the Fig. 2 done for SA Kr (panels a, 
b) and 131 Xe (c,d) targets. 

of parameters and \ 2 in case of symmet- 
ric Lund distribution function are close to 



those for the case of standard Lund distri- 
bution function. Figure 1 shows that the 
good agreement of the leading hadron ap- 
proximation with the data is accidental, be- 
cause the agreement worsens with increas- 
ing of Nh ai i from 1 to the more realistic val- 
ues for multiplicity of hadrons. Compari- 
son with the one dimensional data shows 
that the calculated with the parame- 
ters obtained in result of fit with symmet- 
ric Lund distribution function visually does 
not differ from others. A small difference 
there is in case of two dimensional data 
which are presented in Figs. 2 and 3. Ex- 
perimental points are the two dimensional 
data (filled symbols) from Ref. jl], and the 
one dimensional data (open symbols) from 
Ref. [8|. The theoretical curves are the 
results of fit with the standard Lund dis- 
tribution function - solid curves; with lea- 
ding hadron distribution function - dashed 
curves; with the symmetric Lund distribu- 
tion function - dotted curves. 

4 Conclusions. 

The HERMES data [4] were used to per- 
form the fit for ITSM with symmetric Lund 
distribution function. It was shown that 
using in the fit the distribution of the con- 
stituent formation lengths of hadrons pre- 
sented in the form of finite sum gives the 
stable result beginning with n > 5. For 
the final fit we used the value n — 10. 
Two-parameter's fit gave satisfactory ag- 
reement with data (see Table 1 and Figs. 
2 and 3). Comparison with other versions 
of distribution function used in previous fit 
showed that the results in the cases of the 
symmetric and standard Lund distribution 
functions are close enough. 
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